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The neighborhood of the QCD chiral critical point is characterized by intense fluctuations of 
the chiral field which could, in principle, generate pronounced experimental signatures. However, 
experimental uncertainties which are inherent to heavy-ion collisions, as well as the modest size 
and duration of the formed plasma, will severely attenuate these signatures. Using Monte Carlo 
techniques, we study second-order event-by-event moments of pions as a prototype for signatures 
of the chiral critical point based on the enhancement of the correlation length and event-by-event 
analysis. We test their viability against some realistic ingredients, similar to the ones found in the 
RHIC Beam Energy Scan program. 


I. INTRODUCTION 

The chiral phase diagram of QCD is believed to possess 
the very distinguishing feature of a second-order critical 
end point [2-0| • This point marks the end of a first-order 
transition line and its neighborhood exhibits very unique 
behavior. The discovery of the chiral critical point would 
represent a major breakthrough in the study of the phase 
diagram of strong interactions. Although many mod¬ 
els agree on the existence of a critical point, there is no 
consensus about its location in the (T, he) plane. Fur¬ 
thermore, because of the sign problem, the reliability of 
lattice results is severely compromised for large values of 
the chemical potential. Nevertheless, current estimates 
strongly suggest the chiral critical point might be reach¬ 
able at the current ultrarelativistic heavy-ion collision ex¬ 
periments (HICs), a possibility which brings the hope of 
accessing relevant information about the phase diagram 
of the strong interactions M. 

The existence, location and properties of the chiral 
critical point constitute a crucial issue in the study of 
the phase diagram of strong interactions. Its clarifica¬ 
tion has been the subject of a large amount of theoreti¬ 
cal work and is one of the main goals of the RHIC Beam 
Energy Scan program (BES) ||. 

One candidate source of experimental signatures of this 
point is the increase of long-wavelength fluctuations in 
its neighborhood, and, more specifically, its impact upon 
the event-by-event distribution of observables. Indeed, 
the use of event-by-event correlations of observables as 
signatures of the critical point in HICs was proposed and 
explored in the literature [171 Hot and moments of the 
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proton distribution are already being used in its exper¬ 
imental search Because of critical fluctuations, 

it is expected that these correlations should exhibit non¬ 
monotonic behavior when the freeze-out conditions of the 
plasma formed in the experiments are varied across the 
neighborhood of the critical point, thus providing a pos¬ 
sible signal of its presence. 

However, while in equilibrium a second-order critical 
point is marked by extremely pronounced features, such 
as the divergence of susceptibilities, in realistic condi¬ 
tions equilibrium near this point is hardly attained and 
critical behavior can be dramatically attenuated by both 
dynamical and finite-size effects, especially in the spe¬ 
cific context of the small, short-lived plasma formed in 
HICs fl3 - (l6l[20l - l25| . Moreover, since HICs are very com¬ 
plex experiments, it is not trivial that the nonmonotonic 
behavior emerging from critical fluctuations can be mea¬ 
sured and the viability of such signals should be tested 
in realistic simulations. 

Here, we use Monte Carlo techniques to test second- 
order moments of pions as signatures of the critical point 
in the context of HICs. While estimates of these signa¬ 
tures are available in the literature jTTil il~il | , to our knowl¬ 
edge no previous attempt has been made to test them in 
a more realistic situation using computer simulations. 

In order to introduce some realism, we discuss and in¬ 
clude some sources of spurious contributions which are 
expected in these experiments. Nevertheless, we restrict 
our analysis to a very simple effective theory 0 and 
simplified models for both the spurious fluctuations and 
critical slowing down [23| . Effects from dynamics are con¬ 
sidered only in the estimates of the correlation length and 
finite-size effects are also partially taken into account by 
using Dirichlet boundary conditions on a sphere. We be¬ 
lieve this simplified scenario to be optimistic when com¬ 
pared to the more complex situation of a real heavy-ion 
collision experiment, while taking into account some of 
the most essential features at play in such an environ- 
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ment. 


II. FLUCTUATIONS OF PIONS 
A. Effective theory 

The chiral field a provides an approximate order pa¬ 
rameter for the chiral phase transition of QCD and is, 
for that reason, subject to long-wavelength fluctuations 
in the neighborhood of the critical point. Although this 
field is not directly observable, one would expect these 
fluctuations to have significant effect upon pions, pro¬ 
tons and essentially every field that interacts strongly 
enough with it, increasing event-by-event correlations 
among their corresponding observables 00. „ „ 

We use the same effective theory as Refs. 00 
Since we are interested in long-wavelength fluctuations, 
we adopt a classical approach along with a homogeneous 
approximation for the chiral field: tr(x) = cro- One can 
then define a probability distribution for its zeroth mode, 
cr 0 : 

P[a 0 ] = e ~ Q ^ T , (1) 

where the effective potential f l can be expanded for small 
fluctuations of a 0 (defining (a) = 0): 

al + 0(crl). ( 2 ) 

In Eq. ©, m a is the physical mass of the chiral field, 
which goes to zero at the critical point. Its value includes 
contributions due to interaction with other fields as well 
as thermal effects. 

We are interested in how the fluctuations of cr 0 , given 
by Eq. ©, affect fluctuations of the pions. For this 
purpose, it is appropriate to use the linear sigma model 
and take the coupling between the pions and tro to lowest 
order in the fields, yielding the interaction Lagrangian 

Cint = ~G <j q tt ■ n + O(cj) 4 ), (3) 

which indicates that, in a first approximation, the inter¬ 
action with (Jo has the effect of correcting the pion mass: 

ml = mi 0) 2 + 2 G a 0 . ( 4 ) 

Expression (j4j), together with Eq. ©, defines a Gaus¬ 
sian probability distribution for the pion mass squared, of 
width 2 G £ yjT/V , where the chiral correlation length £ 
is given by m~ 1 . Hence one can calculate event-by-event 
moments of pions by considering an ensemble of realiza¬ 
tions, each with a different pion mass, corresponding to 
chiral field fluctuations. This perspective is especially 
interesting for sampling Monte Carlo events. 

Although Eq. © is not manifestly invariant under chi¬ 
ral transformations, it should be clear that the employed 
coupling is derived from an app roximately invariant po¬ 
tential, as explicit in Ref. [l0| ■ We also note that the 


value of G should depend on medium conditions, being 
estimated at ~ 1900 MeV in vacuum and ~ 300 MeV at 
the critical end point H, although we consider it to be 
constant in the vicinity of the critical point. 


B. Strategy 


In an ultrarelativistic heavy-ion collision with quark- 
gluon plasma formation, this plasma is expected to ex¬ 
pand and cool. At a certain temperature the system 
undergoes chemical freeze-out, and inelastic collisions 
among the particles stop. After further expansion and 
cooling, kinetic freeze-out is achieved, and after some 
particles decay, a resulting distribution of particles is de¬ 
tected. Repeating the experiment many times, one gath¬ 
ers a statistical distribution of observables, which con¬ 
tains information about the thermodynamic properties 
of the system at freeze-out and, hopefully, about whether 
the freeze-out conditions are near criticality. 

We wish to reproduce critical contributions to second- 
order moments of pions in a Monte Carlo simulation. If 
we consider the system to be in equilibrium and the only 
effect of interaction to be a correction Sml to the pion 
mass squared, the state of the pions at freeze-out is com¬ 
pletely specified by the set of occupation numbers {n p } 
where p labels each eigenstate of the free one-particle 
Hamiltonian, including the momentum p and all rele¬ 
vant quantum numbers. Thus, our strategy is to sample 
sets of occupation numbers for each event, corresponding 
to different values of ao or, equivalently, ml , according 
to the effective theory outlined above. This can be ac¬ 
complished by sampling a value for ml from Eqs. © 
and ©, followed by each occupation number, sampled 
independently from the Boltzmann factor e~^^ Up ~^ Up , 
where uj p is the one-particle energy of the mode labeled 
by p and depends on ml. In order to have a finite number 
of modes, it is necessary, of course, to work with a finite 
number of momentum modes, which can be done by im¬ 
posing adequate boundary conditions and introducing a 
cutoff for the momentum. 

It is possible to calculate critical contributions to 
second-order moments within this framework. Because of 
chiral field fluctuations, the correction to an energy level 
uip, Sojp = Sml/2u}p + 0(5ml ), fluctuates simultaneously 
for all modes, introducing correlated fluctuations among 
all occupation numbers. This can be seen by calculating 
the microscopic correlator (A n p A rik) with p ^ fc, where 
An p denotes the fluctuation of n p about its average in a 
given event and (• • •) denotes an average over an infinite 
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number of events, 

(A n p A n k ) = ^ (a n p A n k *V "*') q 

= 4- ((<5"4) 2 > ^ ——— (An p A?z fe n p / n fe /)o+0((<JTr4) 4 ) 
8 p-^ ( w p ' Uk' 

* [ ^F7T--U^ + f P )Ml + h), ( 5 ) 

_/ V Ldp 

which is in agreement with Ref. El and where we have 
used (An p Ank)o = S p k f p ( 1 + f p ), with the subscript o 
indicating nonperturbed averages, Z denoting the parti¬ 
tion function, and f p = — 1) _1 . 

Near the critical point, the correlation length £ grows 
and Eq. (J5]) indicates that second-order moments of pions 
increase quadratically with it, providing possible signa¬ 
tures of criticality. In order to study higher-order cumu- 
lants, couplings of higher order in op should be included 
in Eq. ©. This is done in Refs. [12|, 111 and it is shown 
that these contributions are proportional to higher pow¬ 
ers of £. 


III. PHYSICAL SCENARIO 

The method presented in the previous section allows 
us to sample events that reproduce critical correlations 
among pions. However, this is clearly not enough to 
study the effect of these correlations in a realistic con¬ 
text — useful analysis requires some detail on HICs and 
their relevant underlying background. 

It is hoped that the position of the critical end point 
could be revealed in these experiments through nonmono¬ 
tonic behavior as the freeze-out conditions of the result¬ 
ing plasma, namely its baryonic chemical potential p b 
and its temperature T, are varied around its neighbor¬ 
hood. Chemical freeze-out conditions should be more 
determinant for signatures involving particle multiplic¬ 
ities, while kinetic freeze-out conditions should affect 
signatures related to the transverse momenta spectra. 
Nonetheless, these conditions are not directly controlled 
but rather estimated as functions of parameters such as 
the centrality class of the collision, its center-of-mass en¬ 
ergy yfs and which ions are made to collide. 

Additionally, the nonmonotonic behavior we aim at 
might not be visible in experiments if it consists of a 
very small peak compared to background contributions. 
Its intensity will depend, among other things, on how 
much the chiral correlation length is allowed to grow 
when limited by finite-size and dynamical effects & 
|25| . The relevant background for our analysis comes 
from any source of noncritical fluctuations that affect the 
studied cumulants. This includes the variations of freeze- 
out temperature, chemical potential and plasma volume 
among collisions identified as freezing out in the same 
conditions, which might hide thermodynamical event-by¬ 
event fluctuations. Dynamical effects are only roughly 


taken into account in the estimate of the maximum cor¬ 
relation length reached in the collisions, while finite-size 
effects are only partially included by boundary conditions 
which restrict the possible pionic modes. 


A. Collision parameters 

There are currently no precise predictions for the po¬ 
sition of the critical point in the (T, ps) plane. A recent 
paper [26| claims to have found evidence of the critical 
point at T ~ 165 MeV and p,B ~ 95 MeV, although a 
lower limit to its baryonic chemical potential, p ,e > 450 
MeV, was estimated in Ref. [27] using finite-size scal¬ 
ing. This lower bound is consistent with lattice results 
excluding PbS, 500 MeV and pb/T < 1 for a chiral 
critical point [28|, (29| . Since we want to simulate realistic 
conditions, the idea is to take data from the RHIC Beam 
Energy Scan (BES) and simulate a situation in which the 
critical point is in reach for these experiments. 

Since the considered thermal distribution of pions is 
not directly sensitive to the baryonic chemical poten¬ 
tial ps in our treatment, the importance of experimen¬ 
tal data is to provide a value of the freeze-out tempera¬ 
ture, as well as to allow us to estimate the volume of the 
plasma, as seen in Section HV A I We choose, for that pur¬ 
pose, data from STAR Au + Au collisions at snn = 7.7 
GeV, for which ps = 420 MeV at chemical freeze-out, 
which we consider to be sufficiently high 30 . 

For our simulations, we choose midrapidity and very 
high centrality, focusing at rapidity in the —0.5 < y < 0.5 
range for the 5% most central collisions, with the purpose 
of reducing the effects of anisotropic flow and event-by- 
event fluctuations coming from initial conditions |3ll . I32| . 
Inspired by Refs. [30;, 3j], we use T = 130 MeV for the 
plasma temperature and R p = 6.8 fm for its average ra¬ 
dius in the 0% — 5% centrality class, standing midway 
between chemical (T c h = 145 MeV, R c h = 5.8 fm) and 
kinetic ( Tki n = 116 MeV, I4m = 2000 fm 3 ) freeze-out 
conditions. Since the strongest critical behavior is ex¬ 
pected among the soft pions, and in order to avoid more 
complex behavior in the high transverse momentum re¬ 
gion [3l], [ 32 } , we consider only pions with transverse mo¬ 
mentum pt below 1 GeV. 


B. Critical slowing down 

One of the main limitations for the growth of the cor¬ 
relation length, and hence for the strength of signals of 
criticality, in HICs is the fact that, considering the finite 
duration of the formed plasma, there is not enough time 
for £ to grow by an arbitrarily large factor. Thus, the 
very large increase in the equilibrium correlation length 
while approaching the critical point results in equally 
large equilibration times, necessarily forcing the system 
out of equilibrium, a phenomenon known as critical slow¬ 
ing down [34| . 
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In order to estimate the highest attainable value of the 
correlation length in a collision, we inspire ourselves by 
Ref. [HI, which takes advantage of static and dynamical 
universality class arguments to make rough but robust 
estimates, largely based on qualitative behavior. 

The model of Ref. [H] for the evolution of £ in time t 
can be rewritten in the following way 


d l = A (L\ * 2 3 ~ Z 

dt UoJ U &,(*)/ ' 


( 6 ) 


Eq. © is an educated guess based on the behavior for 
small deviations from equilibrium and depends on the 
universal exponents a = 0.11, v = 0.63 and z = 2 + 
a/v (35L IH|. as well as the dimensionless, nonuniversal 
parameter A, which cannot be directly estimated. The 
model assumes the system to be in equilibrium until it 
reaches, at temperature To, a correlation length £0 large 
enough that universality arguments apply. 

The function £ e<? (i) describes the correlation length for 
a system in equilibrium at temperatures lower than To. 
In order to find it, some simplifying hypotheses must be 
used along with universality. Namely, the plasma is taken 
to cool down at fixed baryonic chemical potential, and its 
trajectory in the phase diagram is supposed to map to 
the Ising model phase diagram so that it is perpendicular 
to its first-order transition line, with the Ising magnetic 
field h being approximated as linear in the temperature 
T. The temperature is also taken to decrease at a con¬ 
stant rate, a simplification which can be improved on, 
although with no significant gain in the results [H]- We 
optimistically consider the system to cool through the 
critical point. From the Ising model universality class, 
£ eq {h) oc where /? = 0.326, <5 = 4.80 [ 2 ! HI, 

and, within these approximations, 


Ceg(f) 


Co 


t 


T 


-v/ftS 


(7) 


where we denote by r := j^° T j the time interval taken 
by the system to cool down from To to Te and choose 
T(—r) = T 0 and T(0) = Te, where Te is the temperature 
at the critical end point. 

Ref- M lacks an estimate for A, a very important 
parameter since it determines how closely £ follows its 
equilibrium value and, consequently, at how high a value 
it peaks. However, it is possible to put an upper bound on 
A by requiring that \dt;/dt\ does not exceed the speed of 
light. The maximum value of \d£/dt\ can be determined 
as a function of A and the ratio x := r/£o- A plot of 
the product x max(|d£/dt|) versus the combination Ax 
is shown in Fig. |T] The resulting curve is very well 
described by x max(|d£/e!t|) = 0.83 (Ax) 0 75 , yielding an 
upper bound A < A max = 1.3 x 0 ’ 33 . 

The peak value of £/£o is controlled by the combination 
A x Q. We choose r = 5.5 fm, inspired by Ref. 33j, 


1 The dependence in the ratio r /£0 can be easily interpreted, since 



FIG. 1. Relation between the unknown constant A and the 
maximum value of \d£/dt\. Restricting the growth of £ to be 
below the speed of light yields the constraint 0.83 (Ax) 0 ’ 75 < 

x. 



FIG. 2. Evolution of the correlation length £ in the most 
favorable scenario, with A = A max . Its value never exceeds 
2 £ 0 , being most likely under 1.8 £ 0 , depending on at which 
instant t freeze-out occurs. 


and, since £0 > 1 /Te, take £0 = 17120 MeV = 1.6 fm, 
yielding x = 3.4 and A max x = 6.60’ □ We consider these 
estimates to be quite optimistic, especially given that our 
value for r is comparable to the lifetime of the system 
[33| . Fig. [2] shows the dependence of £/£o on t/r on the 
most optimistic scenario, namely A = A max . We note 
that, even in this scenario, £ can hardly exceed 1.8 £0 — 
2.9 fm, reaching at most 2.0 £0 = 3.2 fm, its freeze-out 


a larger value of r means a longer time reacting to criticality. 

2 In our case, t — 5.5 fm can be understood as taking Tq — Te = 44 
MeV and \dT/dt\ = 8 MeV/fm, although choosing the value of 
t is less restrictive than fixing these quantities. 

3 The requirement £o > 1 /Te comes from the conditions for the 

application of universality, namely £ 1 /Te, which is never 

exactly satisfied [23ll . 
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value depending on how early freeze-out occurs. 

IV. SPURIOUS FLUCTUATIONS 

Critical behavior is not the only source of event-by¬ 
event correlations in heavy-ion collisions. As the be¬ 
havior one wishes to observe experimentally is marked 
by correlated fluctuations of observables of pions, any 
global fluctuation of experimental parameters is a source 
of background and should hence be taken into account. 
In fact, collisions at the same center-of-mass energy and 
in the same centrality class might correspond to slightly 
different parameters such as freeze-out temperature and 
plasma volume. These fluctuations of the parameters 
have a global effect upon the resulting particles which 
might be indistinguishable from critical collective behav¬ 
ior, affecting the measured correlations and providing 
background to critical fluctuations. 

A. Geometrical fluctuations 

Information about the system size in heavy-ion colli¬ 
sions is usually obtained through centrality binning. This 
implies that events belonging to the same centrality class 
will in general have different volumes, with values within 
a given range. Volume fluctuations are expected to arise 
both from a centrality bin width effect (CBWE), gener¬ 
ated by variations of the volume within a centrality bin, 
and a centrality resolution effect, related to initial vol¬ 
ume fluctuations (37i - l39j . For simplicity, we only consider 
CBWE volume fluctuations. In this case, it is possible 
to estimate the form of the volume distribution corre¬ 
sponding to a given centrality class by considering the 
probability distribution of values of the impact parame¬ 
ter. Furthermore, since the resulting distribution turns 
out not to be Gaussian, it is clear that volume fluctua¬ 
tions will also affect higher-order cumulants. 

Prom a geometrical argument, the number of ways in 
which two nuclei can collide with an impact parameter b 
should be proportional to the perimeter 2tt b of the circle 
defined by it. Therefore, the corresponding probability 
is expected to be linear in b: 

T{b) ocb . ( 8 ) 

Also using simple geometry, and considering the collid¬ 
ing nuclei as discs O it is easy to calculate the transverse 
overlap area A between two equal colliding nuclei as a 
function of b and their radius Rn, 

A(b, R n ) = 2R% cos -1 - b \j R N - 1 -p ( 9 ) 


4 Due to Lorentz contraction, treating the nuclei as discs is more 
reasonable than considering them to be spheres. 


where Rn can be taken from the parameter ro on 
the Woods-Saxon nuclear density profile@ p{r ) oc (1 + 
e (i r- 0 )/a)—i, anc j we uge tire value Rn = To = 6.38 fm 
(from Refs, [iol l4l|~). 

Hence, one can get a probability distribution for the 
volume V of the plasma at kinetic freeze-out by suppos¬ 
ing that this volume is proportional to the initial overlap 
area A, with a proportionality factor C with dimensions 
of length. In order to estimate C , we fix the average value 
of the plasma radius in the 0% — 5% centrality class at 
R p = 6.8 fm, yielding C = 12.7 fm [30l l33|. 

B. Temperature fluctuations 

We know very little about the temperature distribution 
for a given class of collisions and lack a simple model 
that connects geometrical and temperature fluctuations. 
For that reason, we just take a Gaussian distribution of 
temperatures, with a 5% standard deviation. 

Even though the larger heat capacity near the criti¬ 
cal point could diminish the impact of initial conditions 
on the freeze-out temperature, suppressing spurious tem¬ 
perature fluctuations, its growth should also be limited 
by critical slowing down |34l |. Since it should scale as 
C v ~ in equilibrium, where 7 = 1.240 0, [0], we 
just consider that scaling to hold as an approximation 
out of equilibrium and find out that Cy can grow by at 
most 70% — 130% for a very limited interval of time, even 
considering sigma fluctuations to be responsible for 25% 
of the heat capacity when £ = £0 and using £/£o = 2 — 2.5 
& 

C. Additional sources of spurious fluctuations 

Until this point, we have only considered direct pions 
in our analysis. Although pions from resonance decays 
could provide a major source of background for the sig¬ 
natures we study here, our results already place a fairly 
stringent limit on their visibility 0 . Hence, the results 
below should be regarded as the most optimistic scenario 
while taking into account the minimal ingredients of an 
ultrarelativistic heavy-ion collision experiment. Exper¬ 
imental results, on the other hand, could be even less 
compelling. A study including pions from decays is left 
for the future. 

While a more refined analysis would also require tak¬ 
ing longitudinal and transverse flow, as well as the cor¬ 
responding event-by-event fluctuations, into account, the 
implementation of such effects in our Monte Carlo algo¬ 
rithm would imply a large demand for computing power 
and is also left for latter developments. This could be 


5 In Refs. EiTii . [Tfii this density profile is referred to as the two- 
parameter Fermi model (2pF). 
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FIG. 3. Signal as a function of an d £, in proportion to the 
reference value, taken at = 6, £. = 1.6 fm. The variances of 
the charged pion multiplicity, N, and the average transverse 
momentum of charged pions for a single event, Pt, are shown. 
The vertical dashed line, in red, marks the maximum value 
£ = 3.2 fm found in Section fill Bl 



FIG. 4. Behavior of the signal to baseline ratio of the vari¬ 
ance of N as the freeze-out temperature is varied for different 
choices of the parameter A. 


done usin g d ata from blast-wave model fits to experi¬ 
ments [42l - [45l |. 

Since we are concerned with a limited window in phase 
space, we disregard effects from energy conservation. 


V. RESULTS 


chiral correlation length £ at freeze-out. The viability 
of different signatures of criticality could then be tested 
by comparing their values for increasing £ and asking 
whether or not they are sufficiently distinct that some 
kind of nonmonotonic behavior might be detected while 
experimentally probing the critical end point neighbor¬ 
hood. 

Since our model only includes £ through the combina¬ 
tion £ 2 := (G£) 2 and the estimated G and £ are very 
uncertain, we display our results in terms of this quan¬ 
tity. A coupling G = 300 MeV is taken from Ref. p~fSj | . 
where its value near the critical point is extracted from 
considerations using a Ginzburg-Landau effective poten¬ 
tial and the linear sigma model. However, G is expected 
to change significantly according to the medium, growing 
to a value of G ~ 1900 MeV in vacuum M , and it is diffi¬ 
cult to estimate the value of £ 2 far from the critical point. 
For that reason, we use £ 2 = (300 MeV • 1.6 fm) 2 = 6 as 
our baseline value. 

Fig. [3] exhibits our main result, consisting of the rela¬ 
tive deviation from the baseline for the second-order mo¬ 
ments ((A7V) 2 )and ((Apr) 2 ) as functions of £ 2 , where N 
is the charged pion multiplicity and pr is the mean trans¬ 
verse momentum in a given event. Both were normalized 
by the appropriate power of (TV) in each sample in order 
to cancel out system-size dependence. These signatures 
are linear in £ 2 , as should be anticipated from Eq. ©, 
and ((AN) 2 ) is observed to provide the strongest signal 
among them, reachingan increase of almost 10% depend¬ 
ing on the value of £0 The quantity F(pt ), suggested 
as a signature in Ref. 0 was also calculated, showing 
increases of 2.2% for £ = 2.9 fm and 3.0% for £ = 3.2 fm. 

A higher signal, reaching up to 35%, is obtained by 
employing the mixed cumulant (AN Apr). However, 
given the temperature difference of ~ 30 MeV between 
chemical and kinetic freeze-out |3C| , it is very unlikely 
that a significant increase in £ due to criticality will affect 
both the particle multiplicity and transverse momentum 
distributions. Hence, we discard this mixed signature as 
unrealistic. 

By using a linear fit0 it is possible to directly relate the 
signal in ((AN) 2 ) and the value of £ 2 , or, alternatively, 
t/r. Fig. [I]shows how the expected signal depends on the 
difference between the freeze-out temperature T/.o. and 
the critical point temperature Te, assuming Tq — Te = 
44 MeV and using different values for A. For a more 
moderate choice of A, such as A max / 10, the peak in this 
signature is decreased to a mere 1.4%. 

It is unclear whether the tested signatures could be 
visible in an experiment. Although we have found signals 
of almost 8% for the maximum correlation length £ = 
3.2 fm, our choice of parameters lies on the optimistic 
edge of their acceptable ranges — we take the plasma to 


The model, methods and estimates discussed above 
were implemented in a Monte Carlo simulation, and sam¬ 
ples of 10 6 events were analyzed for several values of the 


6 Results were succesfully compared with analytical estimates. 

7 The mentioned linear fit yields good agreement, with a reduced 
X 2 of 2.3. 
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spend almost its entire lifetime near criticality and the 
correlation length to grow almost at the speed of light 
while following its equilibrium value. We remark that 
changing the value of the time scale r has the same effect 
as changing A by the same factor. 

Additionally, albeit it contains some realistic features, 
we have explored a very simplified and optimistic sce¬ 
nario, assuming perfect equilibrium and neglecting fac¬ 
tors such as hydrodynamic flow and its related fluctua¬ 
tions as well as rescattering and resonance decay effects, 
all of which can have a strong impact upon the signal. 
As an equilibrium distribution was considered for pions, 
even though m 2 fluctuates, we have implicitly neglected 
their equilibration times when compared to the time scale 
for chiral critical fluctuations. We take this as a simpli¬ 
fying hypothesis, albeit the equilibration time scales for 
the multiplicity and momenta of pions should further di¬ 
minish the expected signatures, possibly completely ru¬ 
ining the signal in case these time scales are too large. 
One should also notice that fluctuations of the freeze- 
out temperature and baryonic chemical potential should 
blur the sharp peak of Fig. [Jin case it is experimentally 
probed. Bearing all these limitations in mind and using 
Fig. [I]as our reference, we find it rather unlikely that a 
signal, in case there is one, should go above 5% in a real 
experiment. 

To gauge the importance of the spurious signal, the 
results discussed above were repeated without consider¬ 
ing any effects from spurious fluctuations. It was found 
that in this case the variance of the multiplicity N still 
provides the largest signal, although it is almost twice as 
sensitive to the value of £ 2 and reaches almost 15% for 
£ = 3.2 fm. Our results indicate that spurious fluctua¬ 
tions only affect the signatures by dissolving the signal 
and simply add a constant background contribution to 
the signatures, increasing the variance of N, normalized 
by its average, from 1.0 to 2.0 at £ = 1.6 fm. Regardless 
of that, this contribution is larger than the one coming 
from criticality and can definitely make the difference be¬ 
tween a pronounced, detectable signature and one which 
is engulfed by noise. 


VI. FINAL REMARKS 

We have studied second-order event-by-event moments 
of pions as a prototype for signatures of the QCD critical 
end point based on the increase of the correlation length 
in its neighborhood. Those signals depend quadratically 
on the chiral correlation length and are therefore ex¬ 
pected to exhibit nonmonotonic behavior as the neigh¬ 
borhood of the critical point is crossed by experimental 
conditions. In order to test this behavior against the ex¬ 
perimental limitations of HICs, we have made use of sim¬ 
plified models and qualitative arguments, avoiding the in¬ 
troduction of unknown parameters and making estimates 
tending to be moderately optimistic. 

We have found spurious fluctuations to have a signifi¬ 


cant impact upon the tested signatures. The background 
in our simulations, provided by fluctuations of the freeze- 
out conditions, and the limited growth of £, due to crit¬ 
ical slowing down effects, were responsible for estimat¬ 
ing the tested signals of criticality as probably less than 
5%, even though the background was probably underes¬ 
timated and our analysis is not particularly conservative. 

One way to make these signals slightly more pro¬ 
nounced would be to decrease contributions from geomet¬ 
rical fluctuations by using either a more restrictive cen¬ 
trality class or a centrality bin width weighting method 
[38|; |39|. However, in our simple model, geometrical fluc¬ 
tuations contribute only about 30% of the background. 

An especially interesting extension to this work would 
be to test strongly intensive fluctuation measures, such 
as the ^-measure for the transverse momentum [461 1 - as 
signatures of the critical point. These measures are con¬ 
structed from extensive quantities to cancel out not only 
volume dependence but also effects from volume fluc¬ 
tuations and should thus display enhanced performance 
compared to ordinary cumulants [13, HH , although it is 
not clear for us how finite-size effects should affect them. 
However, the separation between kinetic and chemical 
freeze-out conditions makes our methods unreliable for 
combinations of quantities such as multiplicity and to¬ 
tal transverse momentum. A relatively simple solution 
would be to use observables involving different particle 
species [49], but we leave this for future work. 

Our results can also be generalized to second-order mo¬ 
ments of protons, with the difference that the two flavors 
of charged pions would be replaced by the two spin states 
of the proton and quantum statistics would change. How¬ 
ever, while the pion mass is expected to be nearly the 
same at the critical point M , this is not the case for the 
mass of the proton, mp, which is approximately propor¬ 
tional to the chiral condensate and should therefore be 
much smaller near the critical end point [5C| 51]. The 
correction to the mass of the proton due to the chiral 
field fluctuations can be directly extracted from the linear 
sigma model and is 5mp = g cto, yielding 5{rn 2 P ) ~ Gp a o, 
with Gp = gmp [HO, [Hi!] • In the vacuum, / T = 93 MeV 
and m v p c = 938 MeV, resulting in g = m v P ac /f n = 10 
[52j . The rhs of Eq. ([5]) would then change by a factor of 
approximately (Gp/G) (p 2 + m 2 )/(p 2 + m 2 P ), which de¬ 
pends on how g and M P depart from their vacuum values 
near the critical point. It is ~ 1 if g preserves its vacuum 
value and Mp < 30 MeV near the critical point, so that 
protons are not expected to display much stronger signa¬ 
tures. Comparison between sig natures from protons and 
pions can be found in Ref. |l3| . 

Another possible generalization would be to include 
higher-order moments into our analysis, at the cost of 
introducing extra unreliable parameters in our model. 
These moments are expected to display stronger depen¬ 
dence on the correlation length fl2j and experimental re¬ 
sults for their dependence on Wsnn are available in Refs. 
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